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The complete lack of theoretical understanding of the quantum critical states found in the heavy 
fermion metals and the normal states of the high-Tc superconductors is routed in deep fundamen- 
tal problem of condensed matter physics: the infamous minus signs associated with Fermi-Dirac 
statistics render the path integral non-probabilistic and do not allow to establish a connection with 
critical phenomena in classical systems. Using Ceperley's constrained path-integral formalism we 
demonstrate that the workings of scale invariance and Fermi-Dirac statistics can be reconciled. The 
latter is self-consistently translated into a geometrical constraint structure. We prove that this 
"nodal hypersurface" encodes the scales of the Fermi liquid and turns fractal when the system 
becomes quantum critical. To illustrate this we calculate nodal surfaces and electron momentum 
distributions of Feynman backflow wave functions and indeed find that with increasing backflow 
strength the quasiparticle mass gradually increases, to diverge when the nodal structure becomes 
fractal. Such a collapse of a Fermi liquid at a critical point has been observed in the heavy-fermion 
intermetallics in a spectacular fashion. 

PACS numbers: 71.27.-|-a, 71.10.Hf, 71.30.-|-h, 74.72.-h 
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I. INTRODUCTION 

Not long ago, it was taken as physical law that macro- 
scopic systems formed from fermions are Fermi liquids, 
behaving in a scaling sense as the non-interacting Fermi 
gas characterized by the Fermi degeneracy scale. This 
has changed drastically in recent times by the discov- 
ery that in some metals the system of electrons be- 
haves very differently. Clear cut examples are the 'quan- 
tum critical' metallic states found at quantum phase 
transitions in heavy fermion compounds^i^i^i^i^i^, while 
there are indications that the metallic states found in 
the cuprate high-Tc superconductors are of a similar 
kindi^iiiSiiSiiiSiiiii^'i^iii^ By tuning a zero temperature con- 
trol parameter (like pressure, magnetic field, density) one 
encounters Fermi liquids with different Fermi surfaces, 
and the quantum critical regime is found at the zero 
temperature phase transition where one metal turns into 
the other— 1£. Physical properties in this quantum criti- 
cal regime are controlled by powerlaws, indicating that 
the system has become scale invariant, in analogy with 
both thermal phase transitions and the quantum phase 
transitions in bosonic systems. Albeit rooted in quantum 
statistics, the Fermi energy is a scale and has therefore to 
vanish in the quantum critical regime. This is confirmed 
in a spectacular fashion in the heavy-fermion systems: 
the mass of the Landau quasiparticles in the Fermi liq- 
uids on both sides of the transition should be inversely 
proportional to the Fermi energy, and this mass is found 
to diverge to infinity in the quantum critical regime.'^ 

How to think about a fermion liquid without Fermi 
energy? The complete lack of success in understand- 
ing the above phenomena is caused by a deep and gen- 
eral methodological problem in many-particle quantum 
physics. For bosonic problems one can employ the pow- 
erful path integral methods of quantum field theory, di- 
rectly relating e.g. the quantum critical state to the well 
understood statistical physics of classical phase transi- 



tions. For fermionic systems this alley is blocked by the 
infamous minus-sign problem rendering the path integral 
non-probabilistic. In fact the mathematics is as bad as 
it can be: Troyer and Wiesei^ showed recently that the 
sign problem falls in the mathematical complexity class 
" NP hard" , and the Clay Mathematics Institute has put 
one of its 7 one million dollar prizes on the proof that 
such problems cannot be solved in polynomial time. 

Some time ago, Ceperleyii discovered an alternative 
representation for the fermionic path integral that does 
not solve the minus-sign problem in a mathematical sense 
but has as virtue that the fermionic statistics is coded 
in a more manageable way: the 'constrained' fermionic 
path integral. In this framework, the minus signs as- 
sociated with Fermi-Dirac statistics are self-consistently 
translated into a geometrical constraint structure (the 
nodal hypersurface) acting on a residual bosonic dynam- 
ics. Although this nodal surface which contains all the 
data associated with the differences between bosonic and 
fermionic matter is a priori not known for an interact- 
ing fermion problem, in the scaling limit only its aver- 
age and global properties should matter. Henceforth, it 
should be possible in principle to classify all forms of 
fermionic matter in a phenomenological way by classify- 
ing the average geometrical- and topological properties of 
the constraint structure to subsequently use this data as 
an input to solve the resulting bosonic path integral prob- 
lem. This procedure is supposedly a unique extension of 
the Ginzburg-Landau- Wilson paradigm for bosonic mat- 
ter to fermionic matter. Employing the Ceperley path 
integral, in this paper we deliver proof of principle that 
fermion statistics and emergent scale invariance underly- 
ing the critical state can be reconciled. 

This paper is organized as follows. In Section |TT] we in- 
troduce the Ceperley path integral and explain the notion 
of the "nodal hypersurface" . As one anticipates, it has to 
be that the scales of the Fermi liquid are encoded in the 
nodal surface since the residual bosonic dynamics can- 
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not possibly generate these scales by itself. In Section Hill 
we indeed establish a one-to-one correspondence between 
the Fermi degeneracy scale and an average nodal pocket 
dimension and demonstrate that in order for a fermionic 
system to become critical the nodal surface has to turn 
into a scale invariant fractal. This should be regarded as 
our most important result which is further illustrated in 
the remainder of this paper. In Section [IVI we introduce 
the concept of Feynman and Coher^i^ of incorporating 
hydrodynamical backflow effects in a quantum mechani- 
cal wave function and show that fermionic backflow wave 
functions describe many-particle states characterized by 
a hierarchy of increasing number of particle correlations. 
The nodal structures of such wave functions are inves- 
tigated in Section |Vl We find that their nodal surfaces 
change drastically with increasing backflow strength and 
turn into a fractal when the backflow becomes hydrody- 
namical, involving a macroscopic number of particles. A 
detailed fractal analysis is provided in Section IVII In 
Section fVIII we perform Monte-Carlo calculations of the 
momentum distribution to extract the quasiparticle effec- 
tive mass as a function of backflow strength. We find that 
the effective mass diverges exactly at the point where the 
nodal surface turns into a fractal. Finally, in Section lVIIII 
our results are summarized and discussed. 



II. CONSTRAINED PATH INTEGRAL AND 
NODAL SURFACE 

The fermion sign problem becomes apparent when ex- 
pressing the many body density matrix /3i?(R, R'; 
with R = (ri,...,rjv) the position in dN dimensional 
configuration space {d the spacial dimension and N the 
number of particles), and f3 = l/(fcsT) the inverse tem- 
perature, as a path integral over worldlincs {Rr} in imag- 
inary time T (0 < T < hf3), weighted by an action ^[Rt-], 

•J R, — i-T^R,' 

S[Rr] - £'dr{^ill + V{Tlr)}, (1) 

where the sum over all possible A^! particle permutations 
V accounts for the indistinguishability of the particles 
and the alternating sign imposes the Fermi-Dirac statis- 
tics. Here p = par(7') denotes the parity of the per- 
mutation, even permutations enter with a positive, odd 
permutations with a negative sign. The term T^(R) is 
a short hand notation for both external potentials and 
particle interactions. For simplicity, we have considered 
spinless fermions. The partition function is obtained as 
a trace over the diagonal elements of the density matrix, 
Zn{P) = J dRpF(R-, R; corresponding to worldlines 
returning to their starting place or a permutation of it. 
For bosons where the minus signs are absent the partition 
function can be viewed as a classical one describing an 



ensemble of interacting cross-linked ringpolymers. How- 
ever, the fermionic minus signs make this probabilistic 
interpretation impossible. 

Some time ago, Ceperley proved^'' that the fermionic 
density matrix can be calculated as a path integral anal- 
ogous to Eq. II]) but summing only over worldlines that 
do not cross the nodes of the density matrix itself which 
define for each given initial point Rq and inverse tem- 
perature P a {dN — l)-dimensional hypersurface in dN- 
dimensional configuration space, 

17r„,0 := {R|pj.(Ro, R; hf3) = 0}. (2) 

Those hypersufaces act as infinite potential barriers 
allowing only for node avoiding worldines Rr with 
Pf(R, Rr; t) ^ for < r < ;i/3. 

To calculate the partition function we have to integrate 
over the diagonal density matrix pi?(R, R; fi./?) which is 
obtained as a path integral over all worldline configura- 
tion R "> VR which do not cross the nodal surface on 
any time slice and belong to the reach 

r;3(R) = {7 : R ^ R'|pf(R, R(t); t) ^ 0}. (3) 

Because of the anti-symmetry of the fermionic density 
matrix under particle permutations V, 

Pf(R,PR;^/3) = PF{VR,Ii;h(3) 

= {-irpF{R,R;hp), (4) 

all worldline configurations corresponding to odd permu- 
tations have to cross a node an odd number of times and 
are therefore completely removed from the partition func- 
tion. They are exactly cancelled out by all node crossing 
even permutations and we are left with an ensemble of 
all node-avoiding worldline configurations corresponding 
to even permutations, 

pp{K,R;h(3) = — J2 / I?R.e-^[^l/^ (5) 

Remarkably, this representation of an arbitrary 
fermion problem is not suffering from the 'negative prob- 
abilities' of the standard formulation. Surely, one cannot 
negotiate with the NP-hardness of the fermion problem 
and Ceperley's path integral is not solving this problem 
in a mathematical sense. However, the negative signs 
are transformed away at the expense of a structure of 
constraints limiting the Boltzmannian sum over world- 
line configurations. These constraints in turn can be re- 
lated to a geometrical manifold embedded in configura- 
tion space: the 'reach', which is determined by the nodal 
hypersurfaces of the fermion density matrix. This reach 
should be computed self-consistently: it is governed by 
the constrained path integral that needs itself the reach 
to be computed. 
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So far, this self consistent reformulation of the path 
integral has only been used to study Fermi liquids like 
helium-3^^ and many body hydrogen22i2i quantitatively 
within path-integral Monte-Carlo simulations fixing the 
nodal constraints with a free-particle density matrix. At 
a first glance, this might seem to be a very crude approx- 
imation, however, in a scaling sense the Fermi liquid can 
be viewed as a gas of weakly interacting Landau quasi- 
particles. As we will see in the next section, imposing the 
nodal constraint structure of non-interacting fermions on 
the residual bosonic dynamics means nothing but taking 
Landau's Fermi-liquid paradigm for granted. 

Let us further inspect the nodal hypersurface and the 
workings of the constrained path integral. Since the 
fermionic density matrix is odd under particle exchanges 

Pf(Ro,^,jR; m = -pj^(Ro,R; hp), (6) 

the density matrix is zero if two fermions are at the same 
position, Ti = Tj, irrespective of the temperature, the ref- 
erence point Rq, and the interactions between the parti- 
cles. Therefore, the Pauli surface 

P=\J{Il\r, = r,} (7) 

is always a submanifold of the nodal surface r^Ro,/?, 
whereas the dimensionality of the Pauli surface is 
dimP ^ Nd-d. 

In terms of a complete set of fermionic eigenfunctions 
(R) with eigenvalues Ea the fermionic density matrix 
is given by 

PF(Ro,R;n/3) = ^e-'3^°*;(Ro)vI/„(R), (8) 

a 

which, in case of a non-degenerate ground state ^'o(R), 
converges in the limit T —^ to pp(Ro, R; hp — oo) = 
**(Ro)^'(R). Therefore, the nodal surface riR^^^ of the 
finite temperature density matrix becomes independent 
of the reference point Rq in the zero-temperature limit 
and converges to the nodal surface of the ground-state 
wave function, Q {R|vI'o(R) = 0}. 

In Fig. ([1^) a random cut through the nodal hypersur- 
face of the ground-state wavefunction ^'(ri, . . . , tat) = 
A/'det (e*'''""^ )^ ^. of A^ = 49 spinless fermions in 

a two-dimensional periodic box is shown. This partic- 
ular particle number corresponds to a set of momenta 
ki . . . , kjv on a grid Ak — 2tt/L [L the linear dimension 
of the box) forming a closed shell in momentum space and 
therefore to a non-degenerate ground state (see inset of 
Fig. ([1^)). The cut is obtained by fixing TV — 1 parti- 
cles at random positions and tracking down the nodes 
of the wave function moving the remaining particle over 
the system. The algorithms used to find the nodes of 
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FIG. 1: a) Cut through the nodal hypersurface of the ground- 
state wave function of A'' = 49 free, spinless fermions in a 
two-dimensional box with periodic boundary conditions. The 
set of momentum states corresponding to a non-degenerate 
ground state is shown in the inset. The cut is obtained by 
fixing A*' — 1 fermions at random positions (green dots) and 
moving the remaining particle (yellow dot) over the system. 
The nodal surface cut is given by the interface between red 
and blue regions corresponding to negative and positive val- 
ues of the wave function, respectively, whereas absolute values 
are encoded in the color shading. The nodal lines connect the 
A'^ — 1 fixed particles since the [Nd — d)-dimensional Pauli 
surface is a lower dimensional submanifold of the [Nd — 1)- 
dimensional nodal hypersurface. b) Sketch of an allowed 
world-line configuration contributing to the Ceperley path in- 
tegral. On every time slice r one particular particle (yellow 
dot) sees a nodal surface determined by the positions of the 
other particles on this time slice. Due to the meanderings of 
the worldlines of the other — 1 particles (green) the nodal 
surfaces form a 'tent' in space-time attached to the particle 
worldlines. This nodal 'tent' acts as a hard core boundary for 
the worldline of the remaining particle (yellow line). 



free fermion and Feynman backflow wave functions stud- 
ied later on is described in detail in appendix |^ Since 
the Pauli surface is a lower dimensional submanifold of 
the nodal surface the fixed A^ — 1 particles are located on 
the nodal lines. We find that the nodes are very smooth 
forming pockets of the order of the average inter-particle 
spacing r^. An investigation of the nodal structures of 
the finite-temperature density matrissii shows that this 
observation holds at any temperature and that the nodal 
surfaces smoothly approach the ground-state nodes in the 
limit T ^ 0. 

The workings of the constrained path integral is 
sketched in Fig. ([1]). On every time slice r one par- 
ticular particle ri sees the nodal constraint structure 
^^Ro,/3 •= {R|p(Ro, R; t) = 0} determined by the po- 
sitions r2(T), . . . , rjv(T) of the A^ — 1 other particles on 
this time slice. Due to the meanderings of the worldlines 
of the A' — 1 particles the nodal surfaces form a 'tent' in 
space-time attached to the particle worldlines since the 
Pauli surface is always a lower-dimensional nodal hyper- 
surface. Since /9_f(R, R('''); t) ^ for all < t < this 
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nodal 'tent' acts as a hardcore boundary for the world line 
of the particle ri and the particle is not allowed to pene- 
trate or collide with the 'tent'. A world line configuration 
as shown in Fig. ^ does not violate the constraints and 
contributes to the Ceperley path integral. 

In a recent tutorial paper^^ the fermion sign problem 
has been studied within both the conventional signful 
and the Ceperley path integral. It turns out that even 
for the free Fermi gas for which every student in physics 
knows the canonical solution, the constrained path inte- 
gral turns into a highly nontrivial affair. Remarkably, in 
momentum space the constrained path integral directly 
leads to a one-to-one correspondence between the Fermi 
gas and a system of classical atoms forming a Mott in- 
sulating state in the presence of a commensurate optical 
lattice of infinite strength, living in a harmonic potential 
trap of finite strength. This analogy is literal and the 
only oddity is that we are talking about an optical lattice 
system in momentum space. We immediately rediscover 
our canonical picture of the Fermi gas, simply because 
the dynamics of the world lines becomes trivial due to 
the conservation of single-particle momentum. However, 
the workings of the nodal constraints in the real-space 
formulation remain to a great extend puzzling J2 

Let us start with the case of free fermions in one space 
dimension, where the physics of quantum matter can be 
regarded as completely understood.^ The deep reason 
is that quantum statistics has no physical meaning in 
1-l-lD, and it is always possible to find a representation 
where the sign structure drops out completely. It is in- 
structive to find out how this is processed by the Ceperley 
path integral. The special status of the one-dimensional 
case becomes immediately clear since the dimensions of 
the Pauli surface and the nodal hypersurface coincide and 
the two manifold become the same. In this situation it 
becomes quite easy to read the reach. Start out with a 
reference point Rq = {xi^X2, ■ ■ ■ ,xm) ordering the par- 
ticles for instance like xi < X2 < ■ ■ ■ < xm- 'Spread 
out' this configuration in terms of world lines meander- 
ing along the imaginary time direction and the Pauli- 
hypersurface reach tells that only configurations are al- 
lowed where these world lines never cross each other at 
any imaginary time. Therefore, the ordering of the par- 
ticles is preserved on all times and we only have to con- 
sider world-line configurations where every particle re- 
turns to its starting position. The particles become ef- 
fectively distinguishable. A more abstract way of saying 
this is that the Pauli surface (=nodal surface) divides the 
iV-dimensional configuration space into A'^! disconnected 
nodal cells, each corresponding to a certain ordering of 
the particles. We have seen that the Ceperley path inte- 
gral for the one-dimensional Fermi gas is equivalent to the 
one for hard-core bosons in 1-1-1 dimensions or from a sta- 
tistical physics point of view to the problem of an ensem- 
ble of polymers with only steric, hard-core interactions 
in 2 dimensions. This 'Pokrovsky-Talapov' problem^** is 
surely a very serious statistical physics problem since the 
constraints correspond with infinitely strong delta func- 



tion interactions and accordingly everything is about en- 
tropic interactions and order-out-of-disorder physics. Re- 
markably, the one dimensional fermion story can be com- 
pletely understood from this radically statistical physics 
viewpoint^^-^^ by using a 'self-consistent phonon' method 
discovered by HelfricbSZ to deal with the entropic inter- 
actions associated with biological (extrinsic curvature) 
membranes. 

Obviously, in d > 2 where the Pauli surface is a lower 
dimensional submanifold of the nodal hypersurface the 
self-consistency problem inherent to the constrained path 
integral cannot be resolved like in the one-dimensional 
case. This is the reason why a general bosonization pro- 
cedure is lacking in higher dimensions. Although the 
Ceperley path integral has a much richer structure in 
higher dimensions it is surely the case that the higher di- 
mensional Fermi liquids have to know about the 'entropic 
dynamics' characteristic for the (l-l-l)-dimensional case. 
However, a generalization of the Helfrich construction to 
deal with the steric interactions with nodal 'tent' in a 
self-consistent way is lacking. Moreover, in contrast to 
the one-dimensional case, in d > 2 world line configu- 
rations corresponding to even permutations are a-priori 
not ruled out by the constraint structure. Recently, Mi- 
tas demonstrated22i2^ that in the case of free spinless 
fermions the nodal surface has a minimal tiling property: 
it divides the configuration space into two nodal cells 
only. Therefore, it is alway possible to find a continu- 
ous path R T'cvcnR- not encountering a node. In other 
words, all even permutations are on the reach r^(R) and 
contribute to the constrained path integral. 

Phenomenologically the Ceperley path integral can be 
viewed as a bosonic dynamics subject to a geometri- 
cal constraint structure and therefore as a statistical 
physics problem. However, it is a highly non trivial ques- 
tion how to reconstruct the free d-dimensional Fermi gas 
within this framework using the real space representation 
([5]) and only the one-dimensional case can be regarded 
as fully imderstood.^^'^^ From the canonical picture we 
know that the Fermi gas is characterized by a sharp Fermi 
surface at zero temperature and that the thermodynam- 
ics at low temperatures is governed by particle-hole ex- 
citations in the vicinity of the Fermi surface leading for 
instance to a linear specific heat C(T) ^ kBT/Ep for 
T <gC Ep, irrespective of the spatial dimension. This is 
surely unconventional for a system of interacting bosons 
where the spatial dimension enters the exponents of low- 
temperature expansions in a natural way. Recently, it 
was conjectured^^ that the nodal surface constraints act 
in a highly nonlocal way leading to an effective reduc- 
tion of the dimensionality in a way that the systems be- 
have qualitatively like soft-core bosons in 1-l-lD! Such a 
'holographic' principle would also explain why the nodal- 
surface constraints prevent the system from undergoing 
a Bose condensation at finite temperatures and might re- 
late the emergence of a sharp Fermi surface with a con- 
densation exactly at T = 0. 
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III. GEOMETRICAL VIEW ON FERMI-DIRAC 
STATISTICS AND SCALE INVARIANCE 

The constrained path integral is a precise reformula- 
tion of the sign-full path-integral in terms of an effec- 
tive bosonic dynamics subject to a geometrical constraint 
structure and therefore leads the way to a probabilis- 
tic, statistical physics interpretation of fermionic systems 
since the sign structure is absorbed in the nodal hyper- 
surface. A fermionic state like the Fermi liquid is charac- 
terized by scales, the Fermi-energy Ep and momentum 
kp. These scales are alien to any bosonic system. It 
has to be that these scales are uniquely encoded in the 
nodal structure, since the residual bosonic system cannot 
possibly generate these scales by itself. 

In the previous section, we have seen that in dimen- 
sions d> 2 the Pauli surface is a lower dimensional sub- 
manifold of the nodal hypersurface, irrespective of tem- 
perature or the form of interactions. For free fermions 
we find the nodal surface to smoothly connect the lower 
dimensional Pauli surface both at zero temperature (see 
Fig. Il^)) and in the finite temperature case.— There- 
fore, the corresponding nodal structure is characterized 
by a scale, an average nodal spacing of the order of the 
average inter-particle spacing r^. This is clearly seen in 
the two-dimensional nodal surface cut shown in Fig. (JTJi): 
moving one particular particle over the two-dimensional 
box we find nodal lines smoothly connecting the — 1 
other particles, forming pockets with a linear dimension 
of the order of . 

The presence of this scale in the nodal surface can be 
also deduced by a different argument. From the reduced 
one body density matrix which is simply the Fourier 
transform of the single-particle momentum distribution, 
n(r) = Jj^ e*'""nk, given in the limit T ^ by 

n(r) = J dR**(ri,r2,...,rjv)«'(ri-|-r,r2,...,rjv), 

(9) 

one obtains an estimate for the nodal spacing. The ze- 
ros of n{r) correspond to the average displacement of a 
particle to be found on a node. For the Fermi gas the mo- 
mentum distribution in the limit N oo simply turns 
into a step function rik = 6(|k| — kp) with Q{x) = 1 for 
X < and Q{x) = 1 for x > 0. This leads to the Fourier 
transform 



n(r) = Cdikpr) "^'"^ Jd/2{kFr), 



(10) 



where denotes a constant depending on the spatial 
dimension and a Bessel function of the first kind. 
Asymptotically, 

n(r) - (fcj.r)-('^+i)/2 cos (kpr - ^^^r^ , (H) 

signaling a long-range periodicity in the nodal structure 
with an average nodal spacing r„ ~ kp^ ~ r^. This 



asymptotic behavior is also generic for the Fermi liq- 
uid characterized by a discontinuity in nik) at the Fermi 
wave- vector kp- 

From the constrained path integral a one-to-one corre- 
spondence between the existence of a scale in the nodal 
hypersurface and the Fermi energy Ep can be established 
by a simple scaling argument. Let us first assume that 
the nodal hypersurface is characterized by an average 
nodal pocket dimension of the order of the inter-particle 
spacing Ts characteristic for the Fermi liquid. An allowed 
world line configuration is sketched in Fig. ([TJd) where we 
follow the time evolution of the nodal constraint struc- 
ture seen by one particular particle. In the Ceperley path 
integral the world-line configurations R,- are constrained 
by the reach determined by the condition that the world 
lines are not allowed to cross the nodal hypersurface of 
the density matrix itself at any time. From the perspec- 
tive of one particular particle this means that the particle 
has to stay in its nodal pocket at all times < r < ft/?. 
Due to the meanderings of the world lines of the other 
particles the nodal structure seen by the particle fiuc- 
tuates in time leading in a time continuum limit to the 
picture of a nodal surface 'tent' which hangs in space- 
time and acts as a hard-core boundary for the particle. 
Since the Pauli surface is at all times a lower dimensional 
submanifold of the nodal hypersurface the world lines of 
the other particles act as twisted 'tent-sticks' on which 
the 'tent' is hanging. 

Every particle has to stay within the nodal 'tent' which 
is governed by the dynamics of the other particles. At a 
timescale Tc when the average square displacement 



;2(r) = ([r,(r)-r,(0r)=2d— r 

2m 



(12) 



of the worldlines becomes of the order of the average 
nodal spacing, 1{tc) = r„, the particles start to collide 
with the nodal tent. This leads to an average collision 
time Tc = {2d)^^{2m/h)r'^ corresponding to an energy 
scale 



Er 



h 

Tn 



2d ~ 

2m " r2/'^(| -I- 1) 2m 



9/d 



(13) 



where in the last step we have used that the nodal spac- 
ing r„ is of the order of the inter-particle spacing r^, 
r„ ~ Ts and introduced the particle density n — N/V — 
r-'^/Kd with Kd = 7r'^/Vr(d/2 + 1) the volume of the d- 
dimensional unit sphere. We recognize immediately that 
Eq. (jl3[) is just the expression for the Fermi energy in 
d dimensions, Ep E^. In the case of free fermions we 
can even be more quantitative and estimate the average 
nodal pocket dimension r„ as the first zero of the reduced 
one-body density matrix n[r) (|10p which in d dimensions 
is given by r„ = Zd/kp with Zd the first zero of the Bessel 
function Jd/2{z)- This yields Ec = o-dEp with a2 ~ 0.27, 
as ~ 0.30 in two and three dimensions, respectively. 
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Prom the above scaling argument we have learned that 
an average nodal pocket dimension is dynamically re- 
lated to a typical timescale Tc on which the particles 
feel the steric constraints imposed by the nodal surface 
'tent'. The Fermi energy we rediscover immediately as 
the corresponding energy scale Ep = h/Tc- What does 
this imply for a critical fermionic state having no knowl- 
edge whatsoever about Ep as required by the underlying 
scale invariance and as observed in various experiments? 
Turning the above scaling argument around, the absence 
of a Fermi degeneracy scale immediately implies that the 
nodal surface cannot possibly carry a characteristic scale. 
The nodal surface constraints have to act in the same way 
on all time and length scales and therefore, the nodal sur- 
face of any fermionic critical state has to be a scale invari- 
ant fractal. Hence, by using scaling arguments resting on 
the constrained path-integral, we have discovered a phe- 
nomenological principle: The collapse of the Fermi liquid 
at a quantum critical point as observed experimentally for 
instance in the heavy fermion metals is necessarily associ- 
ated with a qualitative change of the nodal surface from a 
smooth to a fractal geometry. This should be regarded as 
the most important finding reported in this Paper since 
it identifies the probabilistic constrained path integral as 
the mathematical framework to reconcile the workings of 
Fermi-Dirac statistics and scale invariance. 



IV. FEYNMAN BACKFLOW WAVE 
FUNCTIONS: PRELIMENARIES 

In the previous section we have convinced ourselves 
that a critical fermionic state is necessarily character- 
ized by an underlying fractal nodal surface. However, 
an inherent difficulty of the present approach is that 
in order to study nodal structures associated with non- 
conventional states a wave-function ansatz is required. 
Let us focus on a concept introduced by Feynman and 
Cohen incorporating hydrodynamical backflow effects in 
a quantum mechanical wave- funct ion 

They argued that the roton in ^He is like a single mo- 
bile atom which is however dressed up by collective mo- 
tions in the liquid. Helium is a nearly incompressible 
fluid in the hydrodynamical sense and the density in the 
neighborhood of the moving particle should be barely al- 
tered. As a consequence there has to be a backflow of 
other particles conserving the total current and leading 
to an enhancement of the effective mass of this quasi- 
particle which can be described quantum mechanically 
by taking a plane-wave wavefunction exp(ikfi) with a 
collective quasiparticle coordinate 

f, = r,+ ^ ?7(r,,)(r, -r,), (14) 

where are the coordinates of the bare particles, rij = 
\ri — Tjl, and r]{r) a smoothly varying function falling 



off like r~^ on large distances corresponding to hydro- 
dynamical, dipolar backflow of an incompressible fluid 
in two dimensions. Much later, it was found out that 
by using fermionic backflow wave functions of the form 

~ det(e'''''"^ ) for node fixing one obtains excellent vari- 
ational energies for the fermionic "^He quantum fluid^ 
and the homogeneous electron ga,sM^ 

It easy to check that such slater determinants of plane- 
wave functions of collective backflow coordinates are in- 
deed obeying Fermi-Dirac statistics since the permuta- 
tion of two particles r^, leads to an interchange of 
the collective coordinates f ^ and rj without changing the 
other collective coordinates and therefore to an overall 
sign change of the determinant. Therefore, as for any 
fermionic state the lower dimensional Pauli surface is 
still a submanifold of the nodal hypersurface of fermionic 
backflow wavefunctions. However, due to the coUective- 
ness build into these wave functions one might expect 
the nodal hypersurfaces to be radically different from 
the smooth free fermion case. This is supported by ear- 
lier work22, reporting a precursor of a roughening of the 
nodal structure in the regime of weak backflow. 

We would like to emphasize that the hydrodynamical 
Feynman backflow build into the fermionic wave func- 
tions has nothing to do with the conventional notion of 
backflow in a Fermi liquid since it involves interactions 
between a macroscopic number of particles at the same 
time. To rationalize this we will derive an expansion of 
the Hamiltonian of which the backflow wave functions are 
eigenstates. In terms of collective backflow coordinates 
we are simply dealing with a gas of free quasiparti- 
cles and the exact Hamiltonian of the system is therefore 
given by 



7i: = ^ek4ck, (15) 

k 

where we have introduced operators cl^, Ck creating and 
annihilating a backflow particle with momentum k, re- 
spectively. To find the representation of this Hamiltonian 
in terms of bare-particle operators cj^, Ck let us first derive 
an expansion of the A'^-particle backflow wave function in 
terms of free-particle states. Using the above Fermi op- 
erators for backflow and bare particles a general relation 
can be written as 



|ki,...,k^)M = cl-.-^ljO) 

= / rqi.---:qivCkj+qj • ■ -Ck^+q^lO) 

Jm,...qN 

' I r 

— / J^qi,...,qjv 
-'qi,...q]V 

x|ki +qi,...,kAr + qAr), (16) 
where the function F is defined as the Fourier transform 
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X ]Je-^'i''''e^^(^o/-('--^), (17) 

i 

and for abbreviation we have defined fi{i^ij) = 
irj{rij)'k.i{Ti — Vj). The general expressions (|16ll7p can 
be further evaluated by expanding the exponentials, 



i hi 

+ ^ /^(l-y)/m(rm„) + ..(18) 



which leads to an expansion of the operator A which 
transforms the free fermion state to the backflow state, 



|ki, . . . ,kjv)bf = A\ki, . . . ,k 



(19) 



Up to second order in the backflow function ij{r) we ob- 
tain 



Ao = 1 



1 



t,j.m,,n 



/»(q)/m(q')cL+q/2Ck 



q "'q' 



^ '^kj Ck^+q'/2'^k„ cj^.^_q,/2Ck„ , (20) 

where we have introduced the Fourier transform 



(q) = 1 1 d'^re-^V.(r) 



(21) 



The Hamiltonian transforms under the inverse of the op- 
erator A, namely 



From the above expansion it becomes immediately 
clear that Feynman backflow generates a hierarchy of pro- 
cesses involving an increasing number of particles. To 
zeroth order in rj, the backflow state is of course triv- 
ially identical to the free fermion one, Aq — 1. Whereas 
the first order terms Ai are solely corresponding to pair 
scattering processes to second order A2 the Hamiltonian 
already contains 3- and 4-body interactions. The ex- 
pansion in powers of 77 (r) is under control only if the 
the strength of the backflow is sufficiently weak. In this 
case backflow can be to a good approximation rational- 
ized in terms of two-body interactions. However, when 
the backflow becomes strong the expansion breaks down 
and n-body interactions of arbitrary order are no longer 
negligible. Surely, such a state characterized by hydrody- 
namical backflow in a true sense involving a macroscopic 
number of particles cannot be obtained perturbatively 
since it requires an inflnite number of vertex corrections. 

Anticipating that we will flnd a critical non-Fermi liq- 
uid state when the backflow turns hydrodynamical, it 
might well be that this gives a way a general wisdom. 
The Fermi liquid is known to be remarkably stable as 
long as one is dealing with a scaling limit where two 
(quasi)particle interactions are less irrelevant than the 
three particle interactions and so forth. The backflow 
ansatz seems to suggest that it is a necessary condition 
for the destruction of the Fermi liquid that the scaling 
flow is such that n-point interactions, with n arbitrarily 
large, are equally marginal. 



V. THE NODAL STRUCTURE OF BACKFLOW 
WAVE FUNCTIONS 

Let us now proceed to calculate the nodal surfaces of 
fermionic Feynman backflow wave functions. In particu- 
lar, we take 



*(ri 



,rw)=AAdet (e^'^'f^: 



(24) 



with a set ki . . . , k^r of momenta corresponding to a non- 
degenerate ground state in the free case (see Fig. (H^)) 
and use collective backflow coordinates ([14]) with a back- 
flow function 



n = (it^ ^noA-\ (22) 

where 'Hq = eiccj^-Ck is the Hamiltonian for the non- 
interacting Fermi gas. To zeroth order (no backflow) we 
trivially obtain A ~ 1 and therefore H = Hq. The next 
corrections to the Hamiltonian in first and second order 
of the backfiow function ri(r) are given by 

Til — ^A^TIq — HqAi 

n2 = AlnoAi + ^I^Aiyno + lnoAl (23) 



^« = (25) 

having the characteristic tail of hydrodynamical 
dipolar backfiow in a two-dimensional incompressible 
fluid. Further, we have introduced a small distance cut- 
off ro and a length a to make the backflow function 77 (r) 
dimensionless and to control the strength of the backflow. 
The expansion in terms of the backflow strength should 
be controlled by the dimensionless parameter a = a/rg 
with the average inter-particle spacing. However, on 
the wave function level there is no restriction to small 
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FIG. 2: Two dimensional cuts through the nodal hypersur- 
faces of fermionic backflow wave-functions for A'' = 49 parti- 
cles and different values of the the backflow strength a = a/rs 
and a small-distance cutoff ro/rs = 0.1. The cuts are obtained 
in the same way as described in the caption of Fig. ([T]). For 
q: = we recover the smooth nodal structure of free fermions. 
With increasing backflow strength, additional clouds of nodal 
pockets start to develop. The linear dimension of these clouds 
scales linearly with a. When the effective backflow range a 
becomes of the order of the inter-particle spacing rs the nodal 
surface qualitatively changes its geometry and seems to turn 
into a fractal. 



values of a and we can in principle follow the evaluation 
of the nodal surface into a regime where the expansion 
breaks down and backflow is governed by infinite number 
of particle correlations. 

Like in the free fermion case we calculate the nodal 
structure on two-dimensional cuts obtained by fixing 
A'^ — 1 particles at random positions and searching for 
nodes when moving the remaining particle over the box. 
Due to the collectiveness of the backflow wave func- 
tions the complexity of the problem increases signifl- 
cantly compared to the free case since the change of the 
coordinate of one particular particle leads to a change of 



all collective coordinates {fi} and we have to recalculate 
the full N X N determinant for every point on the cut. 
For details on the node searching algorithm we refer the 
reader to appendix [Xl 



In Fig. ([2]) the evaluation of the nodal surface on a 
particular cut for = 49 particles with increasing back- 
flow strength a — a/rg is shown. As a reference, in the 
upper left frame the nodal surface cut for free fermions 
corresponding to a = is plotted. As discussed earlier 
the nodal surface smoothly connects the lower dimen- 
sional Pauli surface forming pockets with a characteris- 
tic dimension of the order of the inter-particle spacing 
Vg. For small values of a additional nodal pockets start 
to develop for small particle separations (in the vicinity 
of the Pauli surface) consistent with the local roughen- 
ing of the nodal surface reported previously;^ whereas 
on larger scales the nodal structure looks similar to the 
free case. The size of these clouds of additional nodal 
pockets is much smaller than the inter-particle spacing 
suggesting that backflow is basically governed by two- 
particle correlations as also expected from the expansion 
in terms of powers of ri{r) (see Eq. I16p . This is confirmed 
by the analysis of the nodal structure of two particles 
subject to a mutual backflow. For small particle sep- 
arations we basically observe the same local change of 
the nodal surface as in the A^-particle case. These local 
two-particle effects have surely nothing to do with the 
concept of collective hydrodynamical backflow as origi- 
nally introduced by Feynman.^^ Moreover, a decrease of 
the small-distance cut-off ro leads to the development of 
a local unphysical fractality in the clouds completely un- 
related to collective behavior. A closer inspection of the 
fermionic backflow wave functions used for node flxing in 
Monte-Carlo calculations shows that they belong to the 
regime of weak backflow accounting for local two-particle 
correlations (" exchange-correlation hole" ) in a crude way. 
This is further amplified by the fact that even short range 
backflow functions r]{r) have been commonly usedi^ 



We find that the size of the clouds of additional nodal 
pockets scale linearly with a identifying the parameter 
a as the effective backfiow range (see Fig. ([2])). With 
increasing a, the clouds start to overlap more and more 
forming clusters of interfering backflow patterns signaling 
that the nodal structure can no longer be understood in 
terms of two-particle correlations and that the backflow 
becomes more and more collective. The nodal surface 
seems to develop a scale invariance up to this cluster 
size ^ which can be further amplifled by decreasing the 
small distance cutoff rp. At the point where a becomes 
of order unity the backflow becomes collective involving a 
macroscopic number of particles. At this point ^ becomes 
of the order of the system size and the nodal surface 
seems to turn into a fractal. 
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VI. FRACTAL ANALYSIS 

To demonstrate that we indeed succeeded to produce a 
scale invariant nodal surface, we evaluate the correlation 
integral C (r) which counts the number of pairs of points 
{Xi} on the nodal surface with a separation smaller than 

C{r) = lim ^ 0(|X.-X,|-r), (26) 

n^oo nin — 1) 

i,j — l....n 

where Q{x) denotes the heavyside function, Q{x) = 1 for 
a; < and Q{x) = 1 for a; > 0. For a fractal object the 
correlation integral which is simply the integral of the 
point-to-point correlation function is expected to scale 
as a power law, C(r) ~ r^" with an exponent Dh very 
close to the Hausdorff dimension of the fractal. '^^ 

Of course, it is impossible to map out the full high- 
dimensional nodal surface. Therefore we instead perform 
the fractal analysis on various two-dimensional random 
cuts as shown in Fig. Whereas for these cut pictures 
we have calculated, solely for illustrational purposes, the 
values of the wave function on all points of a fine two- 
dimensional grid, for the fractal analysis we only need 
the points on the nodes which in Fig. ^ correspond to 
the interface between positive and negative regions. To 
track down the nodes on the cut we use a triangulation 
method directly following the nodal lines. This algorithm 
which is described in detail in appendix 1X1 has the virtue 
of omitting the calculation of the highly collective wave 
functions on too many points away from the nodes. To 
test the algorithm we have calculated the nodes on the 
same two-dimensional cut and for the same parameters 
as in Fig. Q and found a perfect agreement between 
the nodal points obtained from the triangulation method 
(see insets of Fig. ([3])) and the interfaces in Fig. 

A fractal nodal surface in full c?A^-dimensional config- 
uration space must have a Hausdorff dimension Dh be- 
tween dN — 1 which is the dimension of a regular nodal 
hypersurface and dN , the dimension of the embedding 
configuration space. Therefore, in d = 2 space dimen- 
sions one expects on a two-dimensional cut a fractal di- 
mension I' = Dh/N in the range 2 — < i/ < 2. 
Hence, with increasing number of particles A'^ we expect 
the Hausdorff dimension v of the nodal surface cut to in- 
crease and to approach v — 2 m. the limit N oo. This 
means that for a large number of particles the nodes on 
the cuts should be very close to space filling. Notice that 
the visual inspection of the cuts can be quite misleading 
regarding the space-filling properties of the high dimen- 
sional nodal surface. 

Let us first analyze the correlation integrals (|26p for 
the nodal surface cuts shown in Fig. To correct for 
finite-size effects we multiply C(r) with a function g{r/L) 
which normalizes the correlation integral by the number 
of available pairs of separation smaller than r in a finite 
box [~-L/2,L/2] X [-L/2,L/2]. This function is given 



a) Oi I b) 




ln{r/rg) ln(r/rs) 



FIG. 3: Correlation integrals C(r) as a function of separation 
in a log-log plot for the same nodal surface cuts and backflow 
strengths a as used in Fig. 0. The corresponding nodal 
surface cuts are shown as insets. The nodal surface turns into 
a fractal at a critical backflow strength Oc ~ 0.8. 

by g{r/L) = 7rr^/(A,.(a;, y))i2, where Ar{x,y) denotes 
the part of the area of a circle with midpoint (x, y) and 
radius r lying within the box. The average {■)l2 has to be 
taken over all points (a;, y) in the box. A straightforward 
calculation yields for r < L/2 

The evolution of the correlation integral is shown in 
Fig. ([3]). Without backfiow {a = 0) where the nodes on 
the two dimensional cut are smooth lines with an aver- 
age spacing of the order of the inter-particle distance 
we find v = I for small distances as expected for a one- 
dimensional object and a crossover around rs to an ex- 
ponent = 2 on larger scales (Fig. ([3^)). This crossover 
signals the existence of an average nodal pocket dimen- 
sion ~ r^: On this scale nodal lines start to see each 
other forming an object which looks two dimensional on 
larger scales. For small backflow strength a we introduce 
another scale a = ar^ in the system which we have iden- 
tified with the size of additional clouds of nodal pockets. 
The change in the small distance behavior can clearly 
be seen in the correlation integral where we find frac- 
tality with a non-universal dimension v up to the scale 
a. On smaller scales, this fractality is cut off by the pa- 
rameter ro <C Ts- At the larger scale Ts we again find 
a crossover to = 2 signaling the existence of an aver- 
age nodal pocket dimension comparable to the free case 
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terns it is possible to calculate C(r) for a large number 
of random two-dimensional cuts. The resulting distribu- 
tions P(y) of the Hausdorff dimensions is shown in the 
inset of Fig. ^ iox N — 13 and = 29 particles, re- 
spectively. For the larger number of particles the Gaus- 
sian distribution becomes narrower and shifts towards 
V = 2. From the mean v and the width of the distri- 
bution Ptq{y) we estimate for the Hausdorff dimension 
of the nodal surface in full diV-dimensional configuration 
space Dh = ND = 2N -1 + Rh with Rh = 0.6±0.3and 
Rh = 0.5 ± 0.3 for = 13 and N = 29, respectively. 



VII. MOMENTUM DISTRIBUTION AND 
EFFECTIVE MASS DIVERGENCE 



FIG. 4: Correlation integrals C(r) for random two- 
dimensional nodal surface cuts for different numbers of parti- 
cles and a critical backflow strength Oc = 0.8. In all cases 
we find scale-invariant behavior and Hausdorff dimensions 
slightly below 2. In the inset we plotted the distributions 
P{v) of the Hausdorff dimensions obtained from a large num- 
ber of random two dimensional cuts for = 13 and N = 29 
particles, respectively. 

(see Fig. (Hb)). Increasing a further, both the range 
of the scale invariant behavior and the fractal dimension 
increase (Fig. ((St))- 

At a critical value ac ~ 0.8 where a becomes com- 
parable to to both scales suddenly disappear and the 
correlation integral C (r) turns into a power law up to the 
system size demonstrating that the nodal surface cut has 
become a scale invariant fractal with a Hausdorff dimen- 
sion v w 1.97 very close to the space-filling dimension 
d — 2 (see Fig. ([SJi))- For values a > ac the nodal 
structure remains fractal without further increase of the 
Hausdorff dimension v. 

We have found that with increasing backflow strength 
a the backflow clouds start to interfere forming clusters 
of linear dimension ^ (see Fig. ([2])). This length scale 
should be identified with the correlation length of scale 
invariant fluctuations corresponding to the crossover in 
C(r) from scale invariant to regular behavior on larger 
scales as indicated by arrows in Figs. (I31d,c). Approach- 
ing the critical value ac, C rapidly becomes of the order of 
the system size and should diverge in the limit N oo. 
However, it is not possible to track down the divergence 
of the correlation length quantitatively since we cannot 
follow the crossover in a regime sufficiently close to ac- 

In Fig. ([4]) we have plotted the correlation integrals 
C{r) for different numbers of particles at a fixed density 
n — N/V ^ rj^ taking the critical backflow strength 
ttc = 0.8. We find power-law behavior over approxi- 
mately three decades for the biggest system and expo- 
nents V slightly below 2 in all cases. Further, we ob- 
serve a small increase of the fractal dimension consistent 
with the inequality 2 — 1/A^ < i/ < 2. For smaller sys- 



We have demonstrated that fermionic Feynman back- 
flow wave functions exhibit very rich nodal structures and 
that by increasing the backflow strength a it is possible 
to continuously tune the nodal surface from the smooth 
one of the free gas to a scale invariant fractal. On the 
other hand, we showed that the Fermi degeneracy scale 
has a simple geometrical meaning: it has a one-to-one 
correspondence to an average nodal pocket dimension. 
This scale disappears as we approach the critical back- 
flow strength. Therefore, we succeeded to produce a crit- 
ical fermionic state of matter lacking a Fermi degeneracy 
scale! 

Experimentally, the collapse of a Fermi-liquid state to- 
wards a quantum critical point has been observed in a 
spectacular fashion in the heavy-fermion metals. From 
Hall resistivity measurements^ it has been concluded that 
the Fermi surface undergoes a discontinuous jump at the 
quantum critical point signaling that there is no sense of 
a Fermi surface underlying the critical state. How does 
the system gets rid of the Fermi energy scale? The quasi- 
particles in the Fermi liquids on both sides of the transi- 
tion are characterized by an effective mass m* ^ 1/Ep, 
which diverges as a power law as one approaches the crit- 
ical point as for example clearly seen in measurements of 
the linear speciflc heat coefficient C{T)/T ~ m*.? A the- 
oretical understanding how the quasiparticles get heavier 
and heavier and finally completely die is completely lack- 
ing and it remains the question to what extend collective 
backflow wave functions can account for such a mysteri- 
ous behavior. 

The quasiparticle effective mass m* is per deflnition 
given by 



with m the bare electron mass and can be expressed by 
use of the Dyson equatio n^^'^^ in terms of the real part 
of the quasiparticle self energy S(k, as 
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m 



_d_ 



ReS(k,w) 



1 



f^ReI](k,^) 



(29) 



k—kp .uj—Q 



Under the conventional assumption that the momentum 
dependence of the real part of the self energy remains 
non-singular a divergence of the quasiparticle mass m* 
can be directly related to a vanishing quasiparticle pole 
strength 



1 - — ReS(k,t^) 



(30) 



k—kp ,0^—0 



which defines the size of the effective Fermi surface dis- 
continuity and the disappearance of the discontinuity Z 
in the quasiparticle momentum distribution ?i(k) would 
correspond to an effective mass divergence m* /m ~ 1/Z . 

Since Z can be derived from the n{k) jump, we calcu- 
late the single-particle momentum distribution of back- 
flow wave functions to establish a connection between 
an effective mass divergence and the disappearance of a 
scale in the nodal hypersurface. From the ground-state 
wave function the single-particle momentum distribution 
at zero temperature can be calculated as a Fourier trans- 
form of the reduced one-body density matrix ([5]), 



n(k) = {2it)-' 
xvl/(ri 



dr / dRe*""**(ri,r2, 



,^n) 



,^n) 



(31) 



For the backflow wave functions this high dimensional 
integral cannot be calculated analytically since the bare 
coordinates enter via the transformations to collective co- 
ordinates p4p in a highly non-trivial way and an expan- 
sion in terms of the backflow strength a (jl6p breaks down 
as we approach the critical value where the backflow be- 
comes highly collective involving a macroscopic number 
of particles at the same time. Of course numerically such 
integrals in D = dN + d dimensions can only be calcu- 
lated by means of Monte Carlo integration. The basic 
idea is very simple and is based on the central limit the- 
orem. Let {xi, . . . ,x„} be a set of uncorrelated points 
in a very high dimensional space which are distributed 
according to a probability distribution P(x) fulfilling the 
requirements -P(x) > and / dxP(x) — 1. Lets intro- 
duce a new random variable Yf = [f{xi) + . . . + f{xn)]/n 
where / is some arbitrary real- valued function with mean 
/i f and variance given by 



A*/ = (/(x))p 

''dx/(x)P(x), 

(/(x) - (/(x)»^ 



J dx(/(x)-M/)'P(x). 



(32) 



(33) 



Then it can be shown that under rather general 
conditional^ that the central limit theorem applies and 
that for large enough n the variable Yf is normally dis- 
tributed with mean fif and standard deviation af/y^, 
irrespective of the dimension. These ideas can easily 
be applied to high dimensional integrals like for the 
momentum distribution we want to calculate. To es- 
tablish a connection with Eq. (|33p we simply rewrite 
the integral by introducing a probability distribution, 
/ = /dx5(x) = /dx/(x)F(x) with /(x) = 5(x)/P(x). 
The integral / can then be estimated by sampling a large 
number n of points Xi according to the probability dis- 
tribution P(x) as 



g(x») 



(34) 



whereas the variance of the estimate of the integral is 
given by 



n 



2 Z-^ 

i=l 



.9(x») 
^(x.) 



n 

n Z — ' 



ff(Xj) 
P(X,) 



(35) 



from which we obtain ±f7// as an estimate for the size 
of the error bar on the computed value of /. Obviously, 
for a given sample size the variance depends significantly 
on the choice of the probability function P(x). In case of 
the integral ri(k) ([?!]) we have tried different distributions 
and found the best convergence for the choice 



P(x) = P(R,r) = l^'(R) 



(36) 



giving for every k point the value of 7i(k) as the sample 
average 



n(k) 



*(ri 



*(ri,r2, . . . ,rAr) 



(37) 



P(R,i 



To sample points R = (ri, . . . , tat) in 2Af-dimensional 
configuration space according to the density profiles 
p(R) = |^'(R)p of backfiow wave functions for differ- 
ent particle numbers and parameters a and rg (see Eq. 
([25)1 ) we use a standard Metropolis rejection algorithm,— 
which has the great advantage that it allows an arbitrar- 
ily complex distribution to be sampled in a straightfor- 
ward way without knowledge of its normalization. For 
a very detailed description of the Metropolis algorithm 
but also for an overview on quantum Monte Carlo meth- 
ods in general we refer the reader to Ref. l40|. Surely the 
density profile p(R) becomes very craggy as we approach 
the critical backflow strength ac ~ 0.8 where the zeros 
of the density form a scale invariant fractal. Since the 
probability distribution of choice (|36p does not depend 
on the coordinates r we just pick these points on the 
plane randomly. 
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FIG. 5: a) Momentum distribution n(k) along kx for differ- 
ent backflow strengths a for A'' = 13 (filled circles), A'' = 49 
(filled triangles), and A'^ = 113 (open squares) particles at a 
fixed density. Here, a relatively large cut-off ro — OArs has 
been used, b) Inverse effective quasiparticle mass m/m* as 
a function of a for ro/r^ — 0.2,0.1,0.05. m/m* has been 
obtained from an interpolation of the discontinuity of n(fc). 
Due to numerical convergence problems we were not able to 
follow the behavior for a > 0.6. Taking the value ac « 0.8 
extracted from the correlation integral analysis the behav- 
ior of n{k) up to a = 0.6 suggests a strong mass divergence 
m* /m '-^ (1 — a/ac)~'' with 7 = 4.5 ± 0.5. 



The resulting normalized momentum distributions 
along the direction for different values of a are shown 
in Fig. ([5^) where the error bars ^Ufj^fn are obtained 
from Eq. ([55)1 . We have used a relatively large small- 
distance cut-off ro — ^Ats to suppress the fractality of 
the nodal surfaces on smallest scales leading to a highly 
oscillatory behavior of the integrands in the ultraviolet. 
For finite size scaling, we have evaluated n(fc) for each 
value of a for various numbers of particles at a fixed 
density leading to an increase the momentum resolution 
Afc = 27r/L. For different numbers of particles we find a 
consistent interpolation of n(k) for all values of a indi- 
cating that the finite size scaling is well behaved. In the 
inset of Fig. ([5^) the corresponding set of momenta en- 
tering the wave functions is shown. To avoid confusion, 
we would like to emphasize that these momentum states 
within the Fermi sphere are fully occupied only in the 
free fermion case. As seen in the expansion (|16p of the 
backflow wave functions in terms of free particle states 
backflow leads to a hierarchy of scattering processes lead- 
ing to a mixing of all kind of excited free particle states 
and must therefore result in a drastic change of the quasi- 
particle momentum distribution. 

For a = we obtain the discrete step function of the 
free Fermi gas. With increasing backflow strength a the 
discontinuity Z decreases continuously consistent with 
the picture that the bare particles get dressed with back- 
flow clouds (see Fig. ([2])) leading to an enhancement of 
the quasiparticle effective mass. In the regime of small 
backflow we find the numerical convergence to be ex- 
tremely fast even for large systems. For larger values of 
a where the backflow clouds start to overlap (Fig. ([2])), 
Z starts to decrease very rapidly. Approaching the crit- 



ical value ttc ~ 0.8 where the nodal hypersurface turns 
into a fractal, for an increasing particle number iV the 
signal to noise ratio goes down significantly. This criti- 
cal slowing down in the numerics indicates that we are 
facing the singular limit of infinite number of particle cor- 
relations. For a > 0.6 it becomes impossible to extract 
the discontinuity since it becomes smaller than the nu- 
merical resolution. Qualitatively, a disappearance of the 
discontinuity at the the critical value ac ~ 0.8 where the 
nodal surface turns into a fractal seems consistent with 
the numerical data. 

To extract the form of the effective mass divergence 
m*/m ~ XjZ as suggested by the disappearing of the 
discontinuity Z we use the value ac = 0.8 extracted form 
the correlation integrals. We find that the behavior of 
XjZ up to a = 0.6 is indeed consistent with a strong 
power-law divergence of the quasiparticle mass 



m 
m 



1 



a 
ac 



(38) 



with an exponent 7 = 4.5 ± 0.5 (see Fig. ([SJd)). Smaller 
values of the cut-off ro lead to a more rapid enhancement 
of the quasiparticle mass for very small values of a but do 
not influence the form of the effective mass divergence. 
This indicates that the effective mass divergence going 
hand in hand with the emergence of scale invariance of 
the nodal structure on large scales is driven by collective 
backflow correlations involving infinite particle interac- 
tions, or in diagrammatic language infinite order vertex 
'corrections'. 



VIII. CONCLUSION 

In summary, employing the constrained path-integral 
formalism we have delivered here proof of principle that 
fermion statistics and the emergent scale invariance un- 
derlying the critical state can be reconciled. We perceive 
it as highly profound that the workings of fermion statis- 
tics in interacting many particle systems can be encoded 
in a geometrical structure (the nodal surface) which in 
turn can be married with the symmetry of scale invari- 
ance to yield a description of fermionic quantum criti- 
cal states. Phenomenologically, the physics of fermion 
systems can be viewed as 'boson dynamics times nodal 
surface geometry' and on this level it is in principle a 
tractable problem to impose a fractal nodal surface, to 
subsequently compute propagators and thermodynami- 
cal properties. Of course, one anticipates that the frac- 
tal dimension of the nodal surface enters thermodynamic 
exponents and sets an anomalous dimension in the prop- 
agators. For the latter we have already found signatures 
from the diffusive behavior imposed by the fractal nodal 
surface generated by backflow. Instead of conventional 
Gaussian diffusion we find super-diffusive behavior corre- 
sponding to single-particle propagators acquiring a Levy- 
flight form. In principle, there is hope for a generalization 
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of Kadanoff-type scaling relations or of Wilsonian renor- 
malization group treatment of fermionic systems around 
such non-Gaussian quantum critical points. 

Also the competition with superconductivity can be 
studied: it has been demonstrated that the nodal struc- 
ture associated with BCS superconductors is subtly dif- 
ferent from that of Fermi gasses^, and it would be quite 
interesting to find out how this affair would work out 
starting from a fermionic critical state. Naively, one ex- 
pects that a fractal nodal surface imposing constraints 
on the dynamics on all length and time scales would lead 
to a drastic enhancement of pairing which might be the 
reason why quite generically instabilities towards super- 
conducting order in the vicinity of fermionic quantum 
critical points have been observed. 

We have used here the backflow wave function just 
as a device to generate a fractal nodal surface, leaving 
open what the actual microscopic conditions are, caus- 
ing criticality in the physical systems. Backflow is as- 
sociated with the effective incompressible nature of the 
fluid flow and this has an interesting resemblance with 
the microscopy of the strongly correlated electron sys- 
tems: in one or the other way, the Mott insulator is close 
by when fermionic quantum criticality is observed and 
'Mottness' renders electron systems to become incom- 
pressible. In fact, the physics associated with backflow 
is quite similar to the interpretation given by Anderson 
of the 'strange metal' as a Gutzwiller projected fermion 
systems, in the sense that both approaches are based on 
a hidden quasiparticle picture and involve singular trans- 
formations of the Fermi liquid state. It needs further 
investigation in what sense the nodal structures associ- 
ated with Gutzwiller-projected wave functions are scale 
invariant and resemble critical states of fermionic matter. 
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APPENDIX A: NODE SEARCHING 

To find the nodes of the free fermion and Feynman 
backflow wave functions we have to track the sign changes 
of ^'(R) on a numerical grid where the nodal struc- 
ture has to look smooth on the scale of the grid size 
e. Obviously, it is not possible to map out the full 
high-dimensional nodal hypersurface in ciA'^-dimensional 
configuration space. Instead, we calculate the nodal 
structures on random two dimensional cuts, obtained 
by keeping N — 1 particles fixed at random positions 
r2, . . . , r^r, and tracking down the zeros of the wave func- 



tion when moving the remaining particle ri over the sys- 
tem [-L/2,L/2] X [-L/2,L/2], 

^^r.,...,r„ = {ri|*(ri,r2,...,rAr) =0}. (Al) 

If the nodal surface is not fractal this defines a 1- 
dimensional hypersurface in the {d — 2)-dimensional con- 
figuration space of the particle ri. To calculate the val- 
ues of backflow wave functions on a certain grid point we 
have to evaluate an (A'' x A^)-dimensional determinant, 
det A with a^j = exp(ikifj) with the collective backflow 
coordinates given in Eq. ([H]) . This is done by using the 
LU-decomposition A — LU where L and U are lower 
and upper triangular matrices of dimension A'^ x A, re- 
spectively. This decomposition is unique if we require 
la = 1 for the diagonal elements of L. In this case, the 
determinant is given by 

det A = det(L-U) = det(L) det(U) = ]J uu. (A2) 

i=l,...,N 

It is easy to check that one can adjust the overall pre- 
factor of the wave function to make it real valued. This 
is a necessary condition to obtain the nodes by the sign 
changes of ^. 

It becomes immediately clear why the collective back- 
flow makes it is much harder to calculate the nodes. In 
the case of free fermions {a = 0) we have = and 
by moving the particle ri over the cut we change only 
the first row of the matrix A. Therefore we can use the 
expansion 

det(A)= J2 (-l)'+'aHdet(AH), (A3) 

i=l,...,N 

where Aii denotes the (A — 1) x (A — 1) matrix obtained 
by dropping the first row and the i-th column of the ma- 
trix A. Only the prefactors an depend on the coordinate 
Ti whereas the sub-determinants have to be calculated 
only once for a particular cut. Including backflow, the 
change of the position ri of the first particle leads to a 
change of all collective coordinates {f i} at the same time 
and the above expansion turns out to be useless. Instead, 
on every grid point we have to recalculate all collective 
coordinates and the full N x N determinant. 

To track the sign changes on the cut we use two differ- 
ent algorithms. The first one is based on a triangulation 
in the vicinity of the nodal lines (see Fig. ([5])) avoid- 
ing the calculation of the wave function on too many 
points away from the nodes. A first point on the node 
is found by following a random walk till we encounter 
a sign change, lets say between the points Ai and Bi. 
We now initialize our triangulation procedure by choos- 
ing one of the two possibilities of completing the side 
AiBi to an equilateral triangle Ai = AiBiCi. If the 
nodal structure is smooth on the scale e of the numer- 
ical grid the nodal line has leave the triangle through 
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FIG. 6: Illustration of the algorithm following a nodal line by 
a triangulation method. 



one of the two new sides BiCi or AiCi, in our example 



it crosses the side BiCi since ^'(i3i)^'(Ci) < whereas 
\l/(Ai)^'(Ci) > 0. We now take the mirror image A2 of 
the point Ai with respect to the side BiCi to obtain the 
next triangle A2 — B1C1A2 and repeat the procedure. 
Using this triangulation we directly follow the nodal line. 

Obviously, the above procedure breaks down when a 
triangle is penetrated by two nodal lines or when a line 
leaves the triangle through the same side it had entered. 
Therefore, the grid size e has to be sufficiently small. In 
particular, even for small backflow the nodal surfaces de- 
velop a local fractality (see Fig. ^) which is suppressed 
on smallest scales by the UV cut-off tq. Hence we have to 
fulfill the requirement e <C to not run into problems. 

For illustrational purposes it is desirable to calculate 
the wave function on all points of a two dimensional nxn 
grid where e — L/n. This method has been used for 
Fig. ^ where the nodal lines correspond to the interface 
between positive (red) and negative (blue) regions. The 
absolute values of 5" are encoded in the color shading. 
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